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ABSTRACT 

A nonlinear  finite  element  analysis  is  performed  for  a finite 
center-cracked  specimen  subjected  to  biaxial  loading.  A Ramberg-Osgood 
type  stress-strain  relation  is  used  to  characterize  the  material 
property.  It  is  found  that  the  energy  release  rate,  J-integral,  stress 
intensity  factor,  strain  intensity  factor  depend  not  only  on  applied 
stress  perpendicular  to  the  crack  but  also  on  applied  stress  parallel 
to  the  crack.  Biaxial  effects  on  fracture  toughness  parameters 
increase  as  applied  stress  increases.  The  coupling  between  biaxial 
effects  and  material  nonlinearity  has  been  indicated. 


1.  Introduction 


Linear  elastic  treatment  of  fracture  is  usually  considered 
applicable  for  net  section  stress  up  to  about  0.8  the  uniaxial 
tensile  yield  stress.  Even  in  this  range  the  immediate  vicin- 
ity of  a crack  tip  has  some  plastic  yield  due  to  the  singu- 
larity of  stress  field,  hence  the  distribution  of  stresses  in 
the  neighborhood  of  a sharp  crack  differs  considerably  from 
what  has  been  predicted  by  linear  elastic  analysis.  Hutchinson 
[1,  2],  Rice  and  Rosengren  [3]  investigated  the  singular  be- 
havior near  the  crack  tip  for  materials  characterized  by  a 
power  hardening  relation  between  stresses  and  strains.  Gold- 
man and  Hutchinson  [4],  Shih  [5,  6]  solved  the  fully  plastic 
crack  problems  and  proposed  some  approximate  functional  rela- 
tions between  J-integral,  stress  and/or  strain  intensity  factor, 
applied  stress,  load  point  displacement,  etc.  Those  works 
just  mentioned  are  based  upon  the  assumption  of  small  scale 
yield.  Hilton  and  Hutchinson  [7]  solved  the  crack  problem 
for  infinite  specimen  by  combining  the  knowledge  of  dominant 
singular  solution  in  the  near-field,  linear  elastic  solution 
in  the  far-field,  and  finite  element  analysis  for  the  region 
in  between. 

As  far  as  the  biaxial  effects  are  concerned,  Lee  and 
Liebowitz  [8]  showed  that,  in  linear  analysis  for  a finite 
center-cracked  specimen,  the  stress  applied  parallel  to  the 
line  crack  does  not  affect  the  stress  intensity  factor  and 
the  biaxial  effect  on  energy  release  rate  is  less  than  0.02%. 

Sih  and  Liebowitz  [9]  showed  that  there  is  no  biaxial  effect 
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on  energy  release  rate  in  the  linear  analysis  of  an  infinite 
center-cracked  specimen.  However,  Kibler  and  Roberts  [10] 
indicated  that  an  increase  in  the  apparent  fracture  toughness 
with  increasing  biaxial  load  was  observed  experimentally,  and 
this  fact  can  not  be  adequately  explained  by  linear  fracture 
mechanics  theory.  This  leads  us  to  believe  the  experimentally 
observed  biaxial  effects  must  be  coupled  with  nonlinearity. 
Hilton  [11]  calculated  the  plastic  stress  and/or  strain  in- 
tensity factor  for  infinite  cracked  plates  subjected  to  bi- 
axial loading.  In  this  work,  we  analyze  a finite  center 
cracked  specimen,  with  Ramberg-Osgood  type  stress-strain  re- 
lation, subjected  to  biaxial  loading,  by  finite  element  method 
to  obtained  fracture  toughness  parameters  such  as  energy  re- 
lease rate,  J-integral,  stress  intensity  factor.  Indeed, 
the  analysis  has  shown  that  the  biaxial  effects  on  those 
fracture  toughness  parameters  are  coupled  with  the  nonlinearity 
introduced  through  the  stress-strain  relations. 


t 
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2.  Stress-Strain  Relations 


In  linear  elasticity,  the  stress-strain  relations  for  a 
homogeneous  isotropic  material  can  be  written  in  either  of 
the  following  forms  [12]: 


°ij  ■ x ekk  5ij  * 2u  Eij  • 


(2.1) 


. . * g‘-  O • • 

i]  2 y 1J 


TT^T)  °kk  6ij  * 


(2.2) 


where  X and  u are  the  Lame  constants,  and  are  the 
stress  tensor  and  strain  tensor  respectively.  Introducing 
stress  deviator  s^  as  follows: 


sij  = °ij  " T °kk  6ij 
(2.2)  can  be  rewritten  as: 


(2.3) 


(1  ♦ *>  °kk*ij  - 


(2.4) 


where  E and  v are  Young's  modulus  and  Poisson's  ratio  respec- 
tively. Adopting  the  model  suggested  by  Ramberg  and  Osgood, 
we  generalize  the  stress-strain  relations  to  the  nonlinear 
range  as  follows: 

Eeij  - (1  ♦ v)sij  ♦ 1 $ 2V  °kk6ij  * 7 aoe  sij»  (2,5) 


where  effective  stress  a is  defined  as: 


2 3 

°e  = 7 sijsij  * 


(2.6) 


For  convenience  sake,  we  introduce  the  following  nondimension - 
alized  quantities: 
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aij 

5 aij/cfY  • 

(2.7) 

= Sij/°Y  . 

(2.8) 

°e 

“ °e/oY  * 

(2.9) 

= Ee^ j / Oy  , 

(2.10) 

n- 1 

v . (2.11) 


Then  (2.S)  can  be  rewritten  as: 


1 2v  3-  n_1- 

£ij  * (1  + v)  i..  ♦ — ^ °kk«ij  + Sij  (2.12) 

In  case  of  simple  tension  test,  namely,  Ojj  * ooy 


and  all  other  ■ 0,  we  obtain  the  following  nonvanishing 
strain  components: 


_n 

■ o ♦ ao  , (2.13) 

1 n 

r22  * F33  " • vo  ' I**  * (2'14) 

The  equations  (2.13)  and  (2.14)  could  be  used  to  determine  the 
numerical  values  of  v,  a,  n experimentally.  A typical  stress- 
strain  curve,  according  to  (2.13),  is  plotted  in  Fig.  1 for 
reference.  In  the  case  of  generalized  plane  stress,  the  stress 
strain  relations  can  be  expressed  in  the  following  matrix  forms 


» 
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3.  Center-Cracked  Specimen 

’•Ve  consider  a rectangular  plate  of  length  2L,  width  2W, 
and  thickness  B,  with  a centered  line  crack  of  crack  size  2a 
subjected  to  symmetric  boundary  conditions  (cf.  Fig.  2). 
Therefore  only  the  first  quadrant  of  the  plate  R = [x,y|0<x<W, 
0<y<L]  needs  to  be  analyzed.  The  boundary  S of  the  first 
quadrant  is  divided  into  five  parts,  i.e.. 
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S -iJ1  S.  , (3.1) 

where  (cf.  Fig.  3) 

S1  - [x-W,  0<y<L]  , (3.2) 

§2  * [y-L,  0<x<W]  , (3.3) 

53  * (x“0,  0<y<L]  , (3.4) 

54  - [y-o,  0<x<a]  , (3.5) 

Sj  ■ [y*0,  a<x<W]  . (3.6) 


We  focus  our  attention  on  the  plane  problem,  in  the  case  of 
generalized  plane  stress,  which  has  the  boundary  conditions 
specified  by: 
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the  biaxial  factor. 


(3.7) 

(3.8) 

(3.9) 

(3.10) 

(3.11) 


The  strain  energy  density  U is  defined  as: 


£ • . 
ij 

U = / a.. 

J 


dEij 


(3.12) 


Utilizing  (2.5  - 2.11),  one  may  obtain 
U = 


1 + v J2  A 1 - 2v  , ,2  an  n + 1 

T-  * —6 <°kk>  * S-TT  °c  /E 
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The  complementary  energy  density  V,  defined  as 


V = 


/ 


£ii  d°i3 


(3.13) 


(3.14) 


is  found  to  be 


UY 

IT 


1 + v — ^ l-2\>  — ^ 

— °e  * 


o -n  + 1 
n + 1 °e 


3 ~ 


(3.15) 


laied  to  eacli  other  as 
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" °ij  eij  - U * (3.16) 

the  total  strain  and  complementary  energies,  per  unit  thickness, 

* * 

U and  V respectively,  of  the  whole  rectangular  plate  have  the 
following  relation 


W 


W L 


'•■/  / V dx  dy  = //  °ij  ui,j  dx  **  - u' 


■W  -L 


-W 


/ 


a.,  n.  u.  ds  - U 
ij  1 i 


(3.17) 

It  is  noticed  that  the  first  term  on  the  right  hand  side  is 
nothing  but  the  line  integral  of  the  inner  product  of  stress 
vector  and  displacement  vector.  Equation  (3.17)  indicates  that 
the  nonlinear  energy  release  rate  G in  fixed  load  and  fixed  grip 
situations  can  be  obtained  as 


9V 

3T7IT 


9U 

9TTay 


Yixed  grip 


(3.18) 


fixed  load 

If  the  following  nondimensionalized  quantities  are  introduced 
as 


x = x/W  , y = y/W  , 


(3.19) 


l = L/W  , c = a/W  , 


(3.20) 
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G can  be  rewritten  as: 


f’« 


dy 


0 0 
1 £ 


fixed  load 


2W 


r1  hj  / u dx  d y 


0 0 


(3.21) 


fixed  grip 

Another  important  quantity  in  fracture  mechanics  is  J- integral 
[13],  It  is  stated  that  the  following  line  integral  J 


J 


(U  dy  - otj  n}  ulx  ds) 


(3.22) 


is  path- independent . Here,  r is  a curve  which  surrounds  the 
crack  tip,  stating  from  the  lower  crack  surface,  traversing 
counterclockwise,  and  ending  on  the  upper  crack  surface,  s 
is  the  arc  length  and  n^  is  outward  unit  vector  norm'l  to  the 
curve.  Because  of  symmetry  and  J being  path  independent,  we 
can  pick  up  any  point  on  S5  as  the  starting  point  and  any 
point  on  S4  as  the  end  point,  and  let  be  any  curve  travers- 
ing counterclockwise  in  R (cf.  Fig.  3),  then  J can  be  re- 
written as: 


2W  a1 


/ 


(tr  dy  - 0tj  n.  Ulj-  ds) 


(3.23) 
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where  the  nondimensionalized  quantities  ui  and  s are  defined 
as : 

oY  W 

= -g — ui  , s = s/W  . (3.24) 

This  problem  will  be  solved  by  using  nonlinear  finite-ele- 
ment method  to  determine  G,  J,  and  other  fracture  toughness 
parameters.  The  procedure  will  be  discussed  in  the  next 
section. 
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4.  The  Procedure 


The  typical  arrangement  of  the  finite  element  layout 
with  212  nodal  points,  i.e.,  424  degrees  of  freedom,  and 
377  triangular  elements  are  illustrated  in  Figs.  4,  5,  6. 

If  j is  the  number  of  a certain  nodal  point,  then  ^2j-l  and 
are  the  nondimensional  displacement  of  the  point  in  x 
and  y direction  respectively,  and  ?2j  are  corre' 

sponding  external  concentrated  force  components  (nondi- 
mensionalized)  acting  on  that  point  (cf.  Fig.  7).  Assume 
the  displacement  field  within  each  element  is  linear  with 
respect  to  coordinates.  This  implies  the  strain,  and 
accordingly  stress,  within  each  element  is  constant.  For 
each  element  let  the  nodal  point  displacements  (I] , strain 
field  [e] , and  stress  field  [a]  be  represented  by 


- 12  - 


then  we  have 

Tel  « [B][T]  , 

[a]  » [D]  [e]  , 


(4.3) 

(4.4) 


where  [D]  is  3 

x 3 matrix  as 

indicated  in 

(2.1 

6) , and 

bi 

0 

bj 

0 

bk 

0 

[B]  ■ 

■k  o 

ci 

0 

C } 

0 

ck 

, (4.5) 

ci 

bi 

c 

j 

bj 

Ck 

bk 

1 

*i 

*i 

2A  - 

det  1 

9 

(4.6) 

1 

xk 

^k 

bi  * 

h • h 

9 

ci  = 

xk 

9 

(4.7) 

with  the  other 

coefficients 

obtained 

by  a 

cyclic  permutation 

of  subscripts  in  the  order  i,  j,  k.  The  stiffness  matrix 
per  unit  thickness  of  this  particular  triangular  element  is 
obtained  as  [14] : 

[k]  - [B]T[D][B]  A , (4.8) 

where  A is  the  area  of  the  element.  The  governing  equation  is 


finally  obtained  as: 
424 


£ Sy  uY  ' • 


B * 1,2,  ....  424 


where  the  424  x 424  stiffness  matrix  is  tha  sum  of  377  local 
stiffness  matrices. 
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Since  matrix  [D]  depends  on  the  effective  stress  ae, 
so  does  matrix  [K] . Thus,  an  iteration  process  has  to  be 
taken  to  solve  the  nonlinear  matrix  equation  (4.9).  And 
also,  it  is  noticed  that  the  principle  of  superposition 
can  not  be  applied.  Therefore,  for  a specifically  given 
applied  stress  a and  biaxial  factor  k,  we  assign  a set  of 
377  trial  values  for  g(I),  I * 1,2,  ...377,  namely,  let 
g(I)  * g*(I)  for  each  triangular  element.  After  solving 
(4.9),  we  have  u*  , then  we  calculate,  for  I = 1,2,  ...377  , 

[e  (I)]  = [B(I)][6(D]  , (4.10) 

[0(1)]  = lD(I)][e(I)]  . (4.11) 

Thus,  the  calculated  values  of  g(I)  are  obtained  as: 

(n-l)/2 

g**(I)  - a[o’(I)  * o*( I)  - ox(I)  ay(I)  + 3o‘y(I)]  (4.12) 

I 

This  iteration  process  will  be  continued  until,  for  each  I, 
the  percentage  difference  between  g*(I)  and  g**(I)  is  below 
certain  allowable  value  of  error.  After  the  iteration  process 
is  completed,  it  is  straightforward  to  calculate  G,  J,  and 
other  quantities  which  we  are  concerned. 
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5.  Crack  Tip  Solution 

In  linear  elastic  fracture  mechanics,  the  stress  inten- 
sity factor  K is  often  taken  as  one  of  the  fracture  criteria 
since  JC  describes  the  singularity  of  the  stress  field  in  the 
neighborhood  of  crack  tip.  For  an  infinite  plate  with  a 
centered  line  crack  subjected  to  uniform  uniaxial  tension, 
the  stress  intensity  factor  and  the  energy  release  rate  are 
related  by  Irwin's  K-G  relation  which  is  represented  as  [9]: 


Kffftc  + 1) 

“"By 


(5.1) 


where  k takes  the  value  (3-4v)  for  case  of  plane  strain  and 
(3-v)/(l+v)  for  case  of  generalized  plane  stress.  In  the 
case  of  nonlinearity  being  introduced  through  stress-strain 
relation  Rice  and  Rosengren  [3],  Hutchinson  [1,2],  Hilton 
and  Hutchinson  [7]  obtained  the  crack  tip  solution  analyti- 
cally. In  this  section  we  recall  some  of  those  results  in 
plane  stress  as  follows  (cf.[l]): 


’ij 


ij 


*. r 


-1 

n+T  .. 


a^ce) 


Y 


-n 

_ n+T  _ 
r eii(6) 


(5.2) 


(5.3) 


n 


where  ® and  the  dimensionless  functions  of  6,  and 

Ejj  , are  detailed  in  [1,2].  In  the  crack  tip  region  the  di- 
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mensionless  J-integral  is  obtained  as  a function  of  the  di- 
mensionless stress  intensity  factor  Y : 


n+1 

a Y c 
o n 


7 , 


(5.4) 


where  cn  depends  on  the  material  hardening  coefficient  n. 

In  plane  stress  cn  takes  the  typical  values  of  3.86,  3.41, 

3.03,  2.87  for  n * 3,5,9,13  respectively.  However,  Hutchinson 
[1],  based  upon  the  assumption  of  small  scale  yielding  and  the 
path  independence  of  J-integral,  set  the  left-hand  side  of 
eqn.(5.4)  equal  to  the  J value  obtained  in  the  linear  and  uni- 
axial case.  In  this  study,  we  do  not  restrict  ourself  in  the 
range  of  small  scale  yielding  and  moreover  we  are  interested 
in  the  biaxial  effects  on  a finite  center-cracked  specimen. 
Therefore,  once  the  finite-element  analysis  has  been  completed, 
we  calculate  the  following  integral, 


y (U  dy  - oi;.  n^  u^-  ds)  = J/2  , 


(5.5) 


along  a curve  traversing  counterclockwise  in  the  first 
quadrant  R and  set 


n+1 

2a  Y_  c 


n ■ / (0  ■ °ij  \ “i.i-  ds) 

IV 


(5.6) 


to  obtain  and  as  a consequence. 
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6.  Results  and  Discussion 


■ 


Given  a specific  center-cracked  specimen,  we  have  a set  of 
material  constants  E,  v,  a,  and  n and  a set  of  geometric  parameters 
of  a finite  rectangular  plate,  namely,  length  2L,  width  2W,  and 
crack  size  2a.  Attention  is  focused  on  case  of  generalized  plane 
stress  with  uniform  biaxial  stresses  o = o,  o = ko  being  applied 

y * 

along  y = +L,  x * +W  respectively.  The  following  nondimensionalized 
quantities  are  introduced: 


* Of  j / Oy  , 

(6.1) 

7ij 

■ Ee^ j / ay  , 

(6.2) 

a = 

n-1 

aOy  , 

(6.3) 

c = 

a / W 

(6.4) 

l = 

L / W 

(6.5) 

X = 

x / W 

(6.6) 

y = 

y / W 

(6.7) 

where  the  yield  stress  oy  is  obtained  according  to  the  usual 

engineering  definition,  namely,  in  simple  tension  test,  when  stress 
o is  equal  to  Oy  , strain  c is  equal  to  oy/E  + 0.002  (cf.  Sechler 


(0.002  E/a) 


1/n 


(6.8) 


Then,  in  the  finite-element  computer  program,  (2.16)  is  used  as  the 
■tress-strain  relation  and  the  boundary  value  problem  is  solved  for 


the  region  R * [x,  y |0  <_  x <_  1,  0<^y<^t]  according  to  the  boundary 
conditions  specified  by  (3.7  - 3.11).  After  the  grand  matrix  equation 
(4.9)  is  solved  by  iteration,  the  following  quantities  will  be  cal- 


culated numerically: 

Complementary  and  strain  energy  densities  (dimensionless)  in  the  I-th 


triangular  element 


J-integral  (dimensionless)  along  any  specified  curve  r 


Strutts  intensity  t actor  (dimension  I ess) 


l/(n+l) 


After  that,  the  nonlinear  energy  release  rate  G,  J-integral  J,  stress 
intensity  factor  K , and  strain  intensity  factor  K are  obtained  as 


G (fixed  load) 


G (fixed  grip) 


1/ (n+1) 


where  A_  is  the  dimensionless  area  of  the  I-th  triangular  element 


For  illustrative  purpose,  we  fix  o*0.02,  n-13,  v-0.33,  lm 2.5 


in  this  work  and  plot  nonlinear  energy  release  rate  G as  function 


of  applied  stress  o for  two  cases,  k«0  and  k**-3,  at  c*0.5  in  Fig. 8 
which  shows  significant  biaxial  effect.  In  other  words,  G in  the 


tension-compression  case  is  higher  than  that  in  the  uniaxial  case 


G,  for  different  values  of  k,  normalized  by  the  linear 


and  uniaxial  energy  release  rate  G,  is  plotted  against  the  applied 


stress.  It  is  noticed  that  for  a/ov  being  less  than  0.3,  the  values 


of  G decrease  as  biaxial  load  factor  k increases,  however,  at  higher 


stress,  one  notices  that  G(k*3)  becomes  even  larger  than  G(k*-1) 


We  believe  this  is  due  to  the  nonlinearity  mainly  caused  by  the  large 


stress  ka  (**1.2a  ) being  applied  along  x*+W.  Also  the  nonlinear  effects 


on  energy  release  rate  are  remarkable.  The  values  of  J-integral  as 


functions  of  biaxial  load  factor  k are  plotted  in  Figs.  10-11  from 


which  we  notice  that  the  general  characteristics  and  the  numerical 
values  of  J-lntegral  are  similar  to  those  of  nonlinear  energy  release 


However,  we  do  detect  that  the  numerical  difference  between  G 
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and  J increases  as  applied  stress  o increases,  especially  at  larger 
k values  (positive  or  negative).  From  Fig. 12  we  see  that  G is  about 
7.8%  lower  than  J when  k=-3  and  it  is  about  4.4%  higher  than  J when 
3 at  a/o  = 0.4  and  c=0.505.  The  biaxial  effect  on  stress  intensity 

factor  is  shown  in  Fig. 13  and, relatively  speaking,  it  is  much  less 
than  that  on  J— integral  and  energy  release  rate  because,  according 
to  eqn.(5.4),  we  have  the  following  relation: 

Ko(k)/Ko(k=0)  = [J<JO/J(k»0)  ]1/(n+1)  . (6.18) 

The  biaxial  effect  on  strain  intensity  factor  is  shown  in  Fig.  14. 
Because  the  dimension  of  linear  stress  intensity  factor  is  different 
than  that  of  nonlinear  stress  intensity  factor,  one  can  not  compare 
these  two  quantities  directly.  Therefore  we  recall  the  definition 
of  small  scale  yield  stress  intensity  factor  (dimensionless) [1] 

kssy  * C<o/oy)2cwAcnl1/(n+1)  • (6.19) 

and  plot  the  percentage  difference  between  K and  K against  o/o„ 

a ssy  3 Y 

for  k * -3,0,3  in  Fig. 15.  Also  we  found  that  even  when  a/oy-*0,  there 

is  still  a 2%  difference  between  K and  K and  this  difference 

a ssy 

is  due  to  the  fact  that  the  linear  value  of  J-integral  for  an  in- 
finite center— cracked  specimen  has  been  used  for  the  right  hand  side 
of  eqn .(5.4). 
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FIG.  4 FINITE  ELEMENT  LAYOUT  (212  NODAL 
POINTS.  377  ELEMENTS) 
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FIG.  6 ARRANGEMENT  OF  FINITE  ELEMENT  LAYOUT  NEAR 
CRACK  TIP  (ENLARGED,  NOT  ON  SCALE) 


Fig.  12  Percentage  Difference  between  Energy  Release 
Rate  and  J— integral 

V*  0.33,  a » 0.02,  n =13,  JL  - 2.5.  O'/ C 1y  ■ 


security  classification  of  this  page  nw>*»  dm  Entered) 


{ REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

1 REPORT  NUMBER 

2.  GOVT  ACCESSION  NO. 

1 RECIPIENT'S  CATALOG  NUMBER 

4.  title  (end  Subtitle) 

THE  NONLINEAR  AND  BIAXIAL  EFFECTS  ON 
ENERGY  RELEASE  RATE,  J- INTEGRAL,  AND 
STRESS  INTENSITY  FACTOR 

*.  TYRE  OF  REPORT  • PERIOD  COVERED 

«.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHOR*,) 

James  D.  Lee  and  Harold  Liebowitz 

I.  contract  or  orant  numbercj 

NAVY  00014-  75-C-  0946 

t.  PERFORMING  organization  name  and  address 

School  of  Engineering  & Applied  Science 
The  George  Washington  University 
Washington,  D.C.  20052 

10.  program  element,  project,  task 

AREA  • WORK  UNIT  NUMBERS 

II.  controlling  office  name  ano  address 

Office  of  Naval  Research 
Arlington,  Virginia  22217 

12.  REPORT  DATE 

March,  19  7 7 

1}.  NUMBER  OF  PAGES 
22 

14.  MONITORING  AGENCY  NAME  • A0DRESSI7/  dllloronl  from  Controlling  Olllce) 

IS.  SECURITY  CLASS,  (of  thle  report ) 

UNCLASSIFIED 

IS*.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

1«.  DISTRIBUTION  STATEMENT  (ol  thle  Report) 

APPROVED  FOR  PUBLIC  RELEASE:  DISTRIBUTION  UNLIMITED 

17.  DISTRIBUTION  STATEMENT  (ol  tho  mbolroel  ontorod  In  »loc*  20.  II  dllloronl  horn  Eoyort) 

t§.  supplementary  notes 

It.  KEY  WORDS  (Contlmto  on  tovotoo  oldo  II  nocoooory  ond  Identity  by  block  numbot) 

Fracture  mechanics  Finite  element  methods 

Biaxial  effects  nonlinear  effects 

Energy  release  rate  J-integral 

Stress  intensity  factor 

\ 

20.  Abstract  (Continue  on  HWN  oldo  II  nocoooory  ond  Identify  by  bleek  number) 

A nonlinear  finite  element  analysis  is  performed  for  a finite 
center-cracked  specimen  subjected  to  biaxial  loading.  A 
Ramberg-Osgood  type  stress-strain  relation  is  used  to  charac- 
terize the  material  property.  It  is  found  that  the  energy 
release  rate,  J-integral,  stress  intensity  factor,  strain 
intensity  factor  depend  not  only  on  applied  stress  perpendicu- 
lar to  the  crack  but  also  on  applied  stress  parallel  to  the  — 

DO  1 » 1473  toirioN  op'i  NOV si  IS  obsolete  UNCLASSIFIED 

S/N  010  2*014*  *401 


SECURITY  CLASSIFICATION  OF  THIS  NAGS  (When  Dote  Entered) 


Ml  T y CLASSIFICATION  of  THIS  PAOEdWnn  P«l»  Enfnd) 


20  ABSTRACT  (continue) 

crack.  Biaxial  effects  on  fracture  toughness  parameters  in- 
crease as  applied  stress  increases.  The  coupling  between  bi- 
axial effects  and  material  nonlinearity  has  been  indicated.^ 

\ 


THE  GEORGE  WASHINGTON  UNIVERSITY 


BENEATH  THIS  PUUUt 
. IS  BUK1KU 

VAULI  FOR  THE  TUTU  RE 
IN  THE  YEAR  O ‘i h 


THT  SIOHY  Ol  INCINH  »IN&  IN  I HIS  HAH  t>!  lilt  I'tAtINt.  l)V  IHE  VAULI 
Five  INI  t RING  HOP  IS  HIM  I HI  MIMUlillOWS  AS  W ft  I I I I.  N IN  lltl.  RECORDS  Ol 
rOUOWIN  r.  r,  OVHINMtNlAL  AND  PRIII  t.SSIllNAL  INCINtUHNO  OHCANliAl  IONS 
' THOSE  Ol  THIS  GEORGE  WASHINGTON  DNIVIRSI1Y 

tll>\UL)  Ol  COMMISSIONERS  1)1  SI  MIC  I Ol  COLUMBIA 
UNITED  SI  VI  IS  ATOMIC  rNTIKjt  COMMISSION 
DlP'RlMFNI  Ol  IHE  ARMS  UNTIED  STATES  Ol  aMIRICA 
DEPART  si  INI  OF  THE  NAVY  liNITFu  SUITS  OF  AMERICA 
» DE  PA  RIMINI  OF  IHE  AIR  TOUCf.  IINl(U)  VI  Alls  Ol  AMERICA 

iK"  national  advisors  coMMirrrv.  tor  afwonauucs 

NATIONAL  BUREAU  OF  STANDARDS  11  S DIPARIMENI  Ol  COMMERCE 
A si  E RICAN  SOC  It  IT  /OF  CIS  I l INCINEIhS 
AMI  RICAN  IN  S T I I t 1 I OF  I 1 I C I R I C A 1.  f N C,  I Ml  RS 
I HI  AMERICAN  SOVIETS  ol  MECHANICAL  INCIMERS 
I HI  SOC  If  r 1 Of  AMERICAN  M I L l l A W N I SCIM  I US 
AMERICAN  INS  I I TE  IE  Ol  MINI  N C.  h M l I A l l E R C.  I C A L I 
0IS1R/CI  Ol  COl.l'  M II I A sorikJV  Ol  PROllSMONAl  E 
nil  INST  I 111  I I Ol  RADIO  IN',  IN  El  RS  INI. 
rill  CHI  *11  CM.  INC  INFERS  CLUB  Ol  WASHINGTON 
W A SH  INCTON  SOCI  t 1 T Ol  I NC  l N | | H S v.  \S 

fail  knir  king  shewn  a s i in  house  - arc  hi  tie  is’' 

CHAR  LFS  H TOMPKINS  COM  t1  ANN  - BUILDERS  „ 

SOClFTY  Ol  WOMEN  ENGINEERS 

' NATIONAL  ACADEMY  Ol  SCIENCES  NATIONAL  R f si.  ARC  H COI  NClIa- 


AND 

1 we 

ARt> 


NC1NII.RS 
•I  O l M F.  R TL.  I HC 


THE  PURPOSE  OF  THIS  VAULT  IS  INSPIRED  BN  AND  DEDTCM  Eb  TO 
CHARLES  HOOK  TOMPKINS.  doct ow  or  t n c,  r s,  r r p v s g* 

fofCATSF  or  HIS  ENGINEERING  CON  1ft  Hi!  LIONS  I <V  T H I S . T Nl  \ F R S IT  N TO  VHS 
wfOMMJ  M.I  YwiO  HIS  NATION  AND  III  ftjfll  !'  T<  \ 1 fh  N ST 

BY  ThE  GEORGE  WASHINGTON  UNIVERSITY. 

tlOHFkl  <?&.ikUlSC,  *"*  •*  ‘"k'*V  f • . Cl  v>Y>  H . M AUVtNt 

— - — sM*  * 


' f*“t  : U *4  A 


wr<?Tp''G- 
■■:T1  - f ■ ’ 


To  cope  with  the  expanding  technology,  our  society  must 
be  assured  of  a continuing  supply  of  rigorously  trained 
and  educated  engineers.  The  School  of  Engineering  and 
Applied  Science  is  completely  committed  to  this  ob- 
jective. 


. 


